Investigation of infrared phonon modes in multiferroic single-crystal FeTe 2 OsBr 
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Reflection and transmission as a function of temperature (5-300 K) have been measured on single 
crystals of the multiferroic compound FeTe20gBr utilizing light spanning the far infrared to the 
visible portions of the electromagnetic spectrum. The complex dielectric function and optical prop- 
erties were obtained via Kramers-Kronig analysis and by fits to a Drude-Lortentz model. Analysis of 
the anisotropic excitation spectra via Drude-Lorentz fitting and lattice dynamical calculations have 
lead to the observation of all 52 IR-active modes predicted in the ac plane and 43 or the 53 modes 
predicted along the b axis of the monoclinic cell. Assignments to groups (clusters) of phonons have 
been made and trends within them are discussed in light of our calculated displacement patterns. 



I. INTRODUCTION 

Recently, significant attention has been given to a 
class of multifcrroics in which specific magnetic ordering 
breaks inversion symmetry, thus allowing spontaneous 
polarization to arise. This magnetically driven ferroelec- 
tric response is the most direct form of magnetoelectric 
coupling and is associated with large magnetoelectrioi 
and magnetocapicative 2 - responses. Electric dipoles in- 
duced by the ordering of magnetic dipoles remain highly 
susceptible to applied magnetic fields, thus making this 
class of multiferroics of interest from a technological point 
of view.— Recent studies on materials such as TbMnO,3, 4 
RbFe(MoC>4)2)^ and CoCrO.^ all show the appearance of 
spontaneous polarization due to non-collinear magnetic 
ordering via spin-orbit mechanisms, viz., the standard 
spin-current and inverse Dzyaloshinski-Moriya models. 
A technological shortcoming of magnetically driven ferro- 
electrics has been the small value of induced spontaneous 
electric polarization observed. Much larger spontaneous 
polarizations have been predicted in materials where fer- 
roclectricity arises from collinear magnetic ordering of 
spins. 1,8 In collinear magnetic ordering, moments may 
vary in amplitude but do not vary in direction; there- 
fore, one does not expect the aforementioned spin-orbit- 
related interactions to be important, 9 enabling the much 
stronger exchange striction interaction as the mechanism 
by which ferroelectricity is induced. 

In this paper, we report our infrared studies on sin- 
gle crystals of the multiferroic FeTe20sBr compound, a 
magnetically driven ferroelectric with nearly-collinear in- 
commensurate spin order. A previous study by Pregelj et 
al*£ has reported significant changes in the magnetic and 
electric orders when external magnetic fields are applied 
along the different crystallographic directions. Strikingly, 
for B || & they observed that fields greater than 4.5 T de- 
stroyed the electric polarization completely. In our inves- 
tigation of magnetoelectric coupling we measured trans- 
mission in the far-infrared down to 15 cm -1 (2 meV) and 
in external magnetic fields up to 10 T oriented along all 



three crystallographic axes. The significant magnetoelec- 
tric coupling previously shown by the application of ex- 
ternal magnetic fields did not surface in the far infrared 
to the extent we examined it. However, we present a 
comprehensive report of the excitation spectrum in single 
crystal FeTe20sBr. The modeling of the infrared active 
phonons as well as lattice dynamical calculations have 
lead us to make an interesting comparison to a previ- 
ous infrared study on a similar geometrically frustrated 
spin-cluster oxyhalide compound, FeTe205Cl~ 1 > where 
an equal number of modes was predicted but a signifi- 
cantly lower number was reported. 

Becker et al*& solved the crystal structure using single- 
crystal x-ray diffraction. The FeTe20sBr compound 
crystallizes in a layered structure where the individual 
layers comprise the be plane in the monoclinic system 
(P2i/c). The layers are weakly connected via van der 
Waals forces and slip at an 18° angle as they stack along 
the normal to the be plane, thus defining the monocline 
a axis. A single layer is composed of [Fe4Oi6] 20 ~ groups 
confined on top and bottom by [Te40ioBr2] 6_ groups. 
The three groups contain common oxygens which connect 
them and also serve to create charge neutrality within 
the layers. The [Fe4Oi6] 20 ~ groups consist of four edge- 
sharing [FeOe] distorted octrahedra. All iron ions possess 
a +3 oxidation state. More detailed descriptions of the 
crystal structure are found elsewherei^ 2 - 

The FeTe205Br compound exhibits nearly-collinear 
(7° of canting between the Fel and Fe2 sites) incommen- 
surate antiferromagnetic ordering of its moments below 
Tjv = 10.6 K i 13 i 14 The amplitude-modulated magnetic 
order is described with the wave vector q=(l/2, 0.463,0). 
Simultaneously, a ferroelectric polarization is induced 
perpendicular to q and the Fe 3+ moments^ The fer- 
roelectric order is attributed to the highly polarizable 
Te 4+ lone pair electrons. Single-crystal x-ray diffraction 
measurements in the ordered state did not detect any 
change in crystal symmetry from the high temperature 
phase; however, clearly distinguishable changes in Fe- 
Te and Fe-0 interatomic distances were observed that 
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signified exchange striction as the means by which the 
inversion center is removed, thus allowing for ferroelec- 
tricity to arise.— Exchange striction is the mechanism 
by which magnetic ions shift away from their centrosym- 
metric positions to maximize their exchange interaction 
energies j^i^ 



II. EXPERIMENTAL PROCEDURES 

Single crystals of FeTe2C>5Br have been grown by stan- 
dard chemical vapor phase method. Mixtures of ana- 
lytical grade purity Fe2C>3, TeC>2 and FeBr3 powder in 
molar ratio 1:6:1 were sealed in the quartz tubes with 
electronic grade HBr as the transport gas for the crys- 
tal growth. The ampoules were then placed horizontally 
into a tubular two-zone furnace and heated very slowly 
by 50°C/h to 500° C. The optimum temperatures at the 
source and deposition zones for the growth of single crys- 
tals have been 500° C and 440° C, respectively. After six 
weeks, many dark yellow, almost orange FeTe205Br crys- 
tals with a maximum size of 8 x 12 x 1 mm 3 were obtained. 
The powder x-ray diffraction pattern obtained by using a 
Rigaku x-ray diffractometer shows the monoclinic space 
group P2i/ c for all FeTe20sBr crystals. 

The zero-field temperature-dependent (5-300 K) re- 
flectance and transmittance measurements were per- 
formed on a 1.3 x 8 x 6 mm 3 single crystal (crystal 1) 
using a Bruker 113v Fourier transform interferometer in 
conjunction with a 4.2 K silicon bolometer detector in the 
spectral range 25-700 cm" 1 and a 300 K DTGS detec- 
tor from 700-7,000 cm -1 . The crystal was cooled using 
a flow cryostat. Room temperature measurements from 
7,000-33,000 cm -1 were obtained with a Zeiss microscope 
photometer. Appropriate polarizers were employed to 
span the desired spectral range. 

Field-dependent transmission measurements in the far 
infrared at a resolution of 0.25 cm -1 were performed on 
a 0.3 x 8 x 5 mm 3 single crystal (crystal 2) at beam 
line U4IR of the National Synchrotron Light Source, 
Brookhaven National Laboratory. The measurements 
employed a Bruker IFS 66-v/S spectrometer in conjunc- 
tion with a 10 T Oxford superconducting magnet and a 
1.8 K silicon bolometer detector. A free standing wire 
grid polarizer was oriented along the preferential polar- 
ization direction of the synchrotron beam. 

At this point it is important to establish with confi- 
dence that we cooled our single crystal below 10 K in 
both the Oxford superconducting magnet and the Janis 
flow cryostat. The superconducting magnet is equipped 
with a variable temperature insert that cools in liquid 
He vapors. This same system has been used to accu- 
rately measure the properties of superconducting mate- 
rials down to 2 Kji£ The accuracy of the cooling capacity 
of the flow cryostat was corroborated by comparing ra- 
tios of T s /Tn for the superconductor Nbo.5Tio.5N with 
those same ratios measured in the superconducting mag- 
net. With the superconductor fastened to the cold finger 



of our flow cryostat with a copper clamp, the cooling 
capabilities, below 10 K, agreed to within ±1 K. 

In materials where the crystal symmetry is lower than 
cubic, electrical excitations strongly depend on the ori- 
entation of the electric field. For these low symmetry 
crystals, the dielectric constant is a 3 x 3 tensor and will 
generally have off diagonal components. For every crys- 
tal symmetry, however, there exists a set of orthogonal 
Cartesian axes, called the principal dielectric axes, which 
diagonalize the dielectric tensorJ^ In an optics experi- 
ment, the principal dielectric axes are found by rotating 
the orientation of the sample which is placed between two 
polarizers that are crossed and remain fixed. Through- 
out a revolution of 360°, four orthogonal orientations will 
result in zero transmitted light beyond the second po- 
larizer. These four dark orientations correspond to the 
electric field of the incoming light aligning itself parallel 
to a principal dielectric axis of the system. 

For a monoclinic system, the crystallographic b axis 
is unique because it is perpendicular to both the a and 
c axes. Accordingly, one principal dielectric axis of the 
system is fixed along the crystallographic b axis, but the 
other two axes can orient themselves arbitrarily as long 
as they remain orthogonal. Consequently, we performed 
single crystal x-ray diffraction on crystal 1 with knowl- 
edge of how the principal dielectric axes fell in relation to 
the crystal facets in order to determine their orientation 
with respect to the a and c crystallographic directions. 

X-ray diffraction measurements were collected at 300 K 
on a four-circle Panalytical X'Pert MRD in a parallel 
beam geometry with Kq radiation. A 9-29 scan was used 
to check that the sample surface corresponded to the be 
plane while pole figure measurements were performed to 
find the spatial orientation of the a axis with respect to 
the be plane. 



III. RESULTS AND ANALYSIS 

A. X-Ray diffraction 

The diffraction experiment showed that the crystallo- 
graphic c axis also coincides with a principal dielectric 
axis of the system. Furthermore, the third principal di- 
electric axis makes an 18° angle with the crystallographic 
a axis. The XRD patterns acquired from the sample with 
the surface aligned perpendicular to the diffraction plane 
are displayed in Fig. [T] Only very narrow (hOO) diffrac- 
tion lines were observed, confirming that the sample was 
a high quality single crystal, with the surface aligned par- 
allel to the be plane. Several acquired pole figures helped 
identify the orientation of the a axis with respect to the 
be plane. The inset to Fig. Q] contains the results of our 
XRD investigation, namely, a sketch of crystal 1 with the 
principal dielectric axes (e^, e^, and ec) superimposed 
on the monoclinic crystal axes. 



3 



4000 



'§ 3000 




30 40 
29 (degree) 

FIG. 1. The XRD pattern for the sample surface aligned per- 
pendicular to the diffraction plane. The inset depicts the ori- 
entation of the principal dielectric axes and crystallographic 
axes with respect to the facets of crystal 1. 



B. Reflectance and Transmittance Spectrum 

The temperature-dependent reflectance spectrum of 
FeTe205Br for electric field oriented parallel Ba, ^b, and 
ec in the frequency interval 30-1,000 cm -1 (4-120 meV) 
is shown in Fig. A strong sharpening of many modes 
accompanied by a hardening of resonance frequencies is 
observed with decreasing temperature. No drastic de- 
viations from the room-temperature spectrum were ob- 
served along any of the three directions when cooling to 
5 K. Particular attention was given to the excitation spec- 
trum slightly above and below the 10.6 K multiferroic 
ordering temperature. The absence of anomalies in the 
spectra upon cooling below 10.6 K gives strong support 
to the observation of no crystal symmetry change in the 
ordered state; however, a subtle contradiction exists^ 
The high temperature P2i/c space group is centrosym- 
metric, but ferroelectricity, for which non centrosymme- 
try is a strict requirement, is observed below T^r. This in 
turn suggests that new infrared modes should accompany 
the transition. However, the small value of polarization 
measured (8.5 /zC/m 2 ) leads one to suspet that a subtle 
transition occurs and additional infrared modes are likely 
below our experimental resolution (<0.5 cm -1 ). 

Figure [3] displays the optical conductivity along the 
three principal dielectric axes in the frequency interval 
30-33,000 cm -1 . The strong anisotropy observed in the 
far infrared persists at high frequencies resulting in three 
distinct absorption edges. The inset of Fig. [3] shows the 
300 K reflectance up to 33,000 cm -1 along Ba, es, and 
ec- 

Crystal 1 exhibits transmission at frequencies above 
the infrared phonon modes and below the absorption 
edge. The dimensions of the crystal only allowed for 
transmission with light polarized along £b and ec- As 
shown in Fig. 21 a sharp dip at 1850 cm -1 as well as a 
broad absorption centered at 3250 cm -1 are present along 
both polarizations. Structures at these particular fre- 
quencies call to mind the absorption bands of H2O, which 
in turn suggest water of hydration in the crystal struc- 
ture. But the anisotropic strengths of these absorptions 



can be used to rule out H2O. Therefore the structures 
suggest absorptions that are intrinsic to FeT^OsBr. The 
differences in energy of the high frequency transmission 
edge along e b and ic reinforce the anisotropy seen in the 
optical conductivity at high frequencies (Fig. [3]). 



C. Field-dependent transmittance 

Crystal 2 transmits at frequencies below the strong 
infrared phonons, in the gaps in-between the phonons, 
and in the mid infrared between the highest phonon and 
the onset of interband transitions. Particular attention 
was given to the transmission below the strong infrared 
phonons. No appreciable change in the spectra were 
observed for transmittance at 5 K measured along €b 
and ec with magnetic fields of 10 T applied along all 
three crystallographic directions. A class of low energy 
magnetic excitations called antifcrromagnetic resonance 
modes (AFMR) are predicted to exist at finite frequen- 
cies in zero field and shift upon the application of exter- 
nal fields. AFMR modes in the FeT^OsBr have been 
reported in a recent electron spin resonance study below 
400 GHz (~13 cm -1 ). 17 However, in an external field of 
10 T the modes remain slightly below our measurable 
frequency range (15 cm -1 ). 



D. Determination of optical properties 

We have used the measured reflectance and transmit- 
tance spectra of crystal 1 to extract the complex dielec- 
tric function. Because Kramers-Kronig analysis requires 
a single-bounce reflectance spectrum over a broad fre- 
quency region, we must first correct the measured re- 
flectance spectrum in regions where the crystal transmits 
by using a combined reflection and transmission analysis. 
An attempt to solve for the single bounce reflectance by 
inverting the full reflection and transmission equations 
leads to a transcendental equation with infinitely many 
solutions. Nevertheless, using the method employed by 
Zibold et o/.,— i.e., assuming that k is small in the re- 
gion of interest, leads to an approximation of reflection 
and transmission that can be solved numerically for a 
single solution. 

The measured reflectance (not shown) and transmit- 
tance (Fig. 2]) in the frequency interval 1500-10000 cm -1 
were significantly lower than expected throughout this 
region of low absorption (Fig. [3]). The low measured val- 
ues were attributed to the scattering of light from facets 
existing on the back surface of crystal 1. To compen- 
sate, we assumed a mask on the back surface of crys- 
tal 1 which lead to scaling the measured reflectance and 
transmittance up by a factor of 1.4 in the interval 1500- 
10000 cm -1 , thus explaining the discrepancy between 
Fig. |4] and Fig. [5j The scaled reflectance and transmit- 
tance yielded a single bounce reflectance which agreed 
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FIG. 3. (Color online) The optical conductivity in the fre- 
quency interval 30-33,000 cm - . The inset shows the 300 
K broadband reflectance out to 33,000 cm -1 from which the 
optical conductivity was extracted via Kramers-Kronig rela- 
tions. 



well with the measured bulk reflectance in regions where 
crystal 1 did not transmit, as shown in Fig. [5] 

Kramers-Kronig analysis can be used to estimate 
the real and imaginary parts of the dielectric function 
from the bulk reflectance R(uj)J£ Before calculating the 
Kramers-Kronig integral, the low frequency data (0.1— 
30 cm -1 ) were extrapolated using parameters from the 
complex dielectric function model described in Sec. Ill E. 
At high frequencies (8 x 10 4 -2.5 x 10 8 cm -1 ) the re- 
flectance was approximated from known x-ray scatter- 
ing functions of the constituent atoms. The gap be- 
tween our highest measurable frequency and the x-ray 
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FIG. 4. (Color online) Transmittance of crystal 1 in the mid- 
infrared and near-infrared regions. 



reflectance data was bridged with a Lorentz oscillator at 
60,000 cm -1 . Subsequently, the phase shift on reflection 
was obtained via Kramers-Kronig analysis, and the op- 
tical properties were calculated from the reflectance and 
phase. 



E. Oscillator-model fits 

The measured reflectance was fit with a Drude-Lorentz 
model to obtain a second estimate of the complex dielec- 
tric function in the infrared range. The model assigns 
a lorentzian oscillator to each distinguishable phonon in 
the spectrum plus a high frequency permittivity, e^, to 
address the contribution of electronic excitations. The 
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FIG. 5. (Color online) Single-bounce reflectance (red) as ap- 
proximated from the scaled measured reflectance (black) and 
the scaled measured transmittance (blue) of crystal 1. The 
elevated level of measured reflectance throughout the trans- 
mission interval is expected due to additional light reflected 
from the back surface of the crystal. 



model has the following mathematical form: 



3=1 



Srf 



(1) 



where Sj , ujj , and jj signify the oscillator strength, center 
frequency, and the full width at half max (FWHM) of 
the jth lorentzian oscillator. The Drude-Lorentz complex 
dielectric function is used to calculate the reflectivity, 
which is compared to the original measured quantity in 
Fig. [6j Similar qualities of fits were obtained for all other 
temperatures and polarizations. 



TABLE I. Oscillator parameters for experimentally observed 
phonons along e.A (at 20 K). The motions of atoms associated 
with groups of phonons have been assigned when the pattern 
was apparent. The atoms are listed in descending order of 
their calculated net atomic displacement. 



Idx 


Osc Str 


Ctr Freq 


FWHM 


Degen 


Asgn 




S 


uj (cm -1 ) 


7 (cm -1 ) 


idx in ec 




1 


0.207 


45.5 


2.0 


1 


Br 


2 


1.320 


79.6 


2.1 


5 


- 


3 


0.356 


88.4 


1.1 


6 


Br,Te,0,Fe 


4 


0.042 


98.1 


1.0 


8 


Br,Te,0,Fe 


5 


0.056 


107.5 


2.1 


9 


Br,Te,0,Fe 


6 


0.048 


137.1 


1.9 


12 


Br,Te,0,Fe 


7 


0.015 


157.6 


1.6 


14 


Te,Br,0,Fe 


8 


0.023 


166.0 


2.2 


15 


Te,Br,0,Fe 


9 


0.354 


184.9 


1.3 


16 


Te,0,Fe 


10 


0.519 


194.3 


2.4 


17 


Te,0,Fe 


11 


0.159 


211.7 


2.0 




Te,0,Fe 


12 


0.007 


236.4 


4.9 


22 


- 


13 


0.084 


247.8 


1.7 


23 


- 


14 


0.009 


256.5 


4.0 


24 


- 


15 


0.169 


291.9 


1.8 


27 


Fe,0,Te 


16 


0.030 


302.8 


3.4 




Fe,0,Te 


17 
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321.6 


3.6 


28 


Fe,0,Te 


18 
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19 
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2.0 


29 


Fe,0,Te 


20 


0.102 
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0,Fe,Te 
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35 


0,Fe 
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32 
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FIG. 6. (Color online) Experimental reflectance (blue points) 
and calculated reflectance (red line) from the Drude-Lorentz 
model of the FeTe20sBr 100 K e_g dielectric function. 



IV. DISCUSSION 

A. Group theory and lattice dynamics 

Due to the lack of anomalies upon cooling below 
10.6 K, we turn our focus to the highly anisotropic far- 
infrared excitation spectrum. The number of phonon 
modes expected in FeT^OsBr can be found by group- 
theory analysis. Using the SMODES program^ we ar- 
rive at the following distribution of modes: 

T o P ucai _ MA (R) + MB (R) + 53^) + ; (2) 

where (R) and (IR) denote Raman active and infrared ac- 
tive modes respectively. The 53A U modes are expected 
to lie along the unique b axis (e^) of the monoclinic sys- 
tem. The 52B U modes are expected to lie in the ac 
plane, and one should be able to resolve all 52 modes 
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TABLE II. Oscillator parameters for experimentally observed TABLE III. Oscillator parameters for experimentally ob- 

phonons along ec (at 20 K). The motions of atoms associated served phonons along e_e (at 20 K). The motions of atoms 

with groups of phonons have been assigned when the pattern associated with groups of phonons have been assigned when 

was apparent. The atoms are listed in descending order of the pattern was apparent. The atoms are listed in descending 

their calculated net atomic displacement. order of their calculated net atomic displacement. 
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7 (cm -1 ) 


idx in &a 






S 


ui (cm - ) 
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1 


4.275 


46.7 


0.2 


1 


Br 


1 


1.388 


35.7 


0.8 


Br 


2 


1.335 


53.9 


0.5 




Br 


2 


0.330 


43.5 


1.4 


Br 


3 


0.951 


65.8 


1.0 




- 


3 


0.899 


53.2 


0.8 


- 


4 


0.216 


71.7 


0.8 




- 


4 


0.289 


60.5 


0.9 


Br,Te,0,Fe 


5 


0.435 


81.2 


0.9 


2 


- 


5 


0.119 


79.1 


1.3 


Br,Te,0,Fe 


6 


0.205 


89.3 


0.9 


3 


Br,Te,0,Fe 


6 


0.075 


85.0 


1.2 


Br,Te,0,Fe 


7 


0.061 


97.0 


0.8 




Br,Te,0,Fe 
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0.052 


91.2 


1.0 


Br,Te,0,Fe 
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1.1 


4 


Br,Te,0,Fe 


8 
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93.3 


1.3 


- 
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0.533 
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0.8 


5 


Br,Te,0,Fe 
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0.9 


- 


10 


0.294 


112.8 


1.0 




Br,Te,0,Fe 
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1.359 


102.0 


1.1 


- 


11 


0.369 


128.5 


0.8 




Br,Te,0,Fe 


11 


0.073 


106.9 


0.7 


- 


12 


0.127 


135.4 


0.9 


6 


Br,Te,0,Fe 
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0.170 


123.7 


0.7 


Te,Br,0,Fe 


13 


0.411 


140.6 


0.9 




- 


13 


0.092 


126.3 


0.9 


Te,Br,0,Fe 


14 


0.857 


157.0 


0.5 


7 


Te,Br,0,Fe 


14 


0.074 


131.4 


1.1 


Te,Br,0,Fe 


15 
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166.8 


1.0 


8 


Te,Br,0,Fe 


15 


1.066 


184.0 


0.4 


Te,0,Fe 


16 


0.212 


184.5 


0.6 


9 


Te,0,Fe 


16 


0.128 


190.2 


0.7 


Te,0,Fe 


17 


0.039 


194.8 


1.2 


10 


Te,0,Fe 
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0.491 
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0.8 


Te,0,Fe 
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204.0 


1.4 
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18 
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0.3 


- 


19 
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19 
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215.3 


1.3 
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20 


0.085 


221.7 


1.1 




- 


20 


0.053 


222.0 


1.1 


- 


21 


0.101 


225.4 


0.9 




- 


21 


0.058 


240.1 


1.3 


- 


22 


0.065 
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1.0 


12 


- 


22 


0.028 


249.0 


0.9 


- 


23 


0.107 


248.0 


1.5 


13 


- 


23 


0.274 


273.9 


0.9 


Fe,0,Te 


24 


0.291 


254.2 


0.6 


14 


- 


24 


0.102 


275.8 


1.4 


Fe,0,Te 


25 


0.239 


274.4 


1.1 




- 


25 


0.661 
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1.3 


Fe,0,Te 


26 


0.195 


284.9 


1.2 




- 


26 
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1.1 


Fe,0,Te 


27 


0.633 


294.8 


1.4 


15 


Fe,0,Te 


27 


0.101 


324.3 


1.5 


Fe,0,Te 


28 
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1.9 


17 


Fe,0,Te 


28 


0.261 


341.7 


1.2 


Fe,0,Te 


29 
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2.2 


19 


Fe,0,Te 
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0.808 
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2.6 


Fe,0,Te 


30 
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426.7 


2.4 


24 


- 


30 
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418.0 


4.1 


- 


31 


1.562 
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using the two orthogonal polarization spectra measured 
in this plane. However, the merging of weaker modes 
near stronger modes can hinder this count. We therefore 
employed lattice dynamical calculations to determine the 
relative resonance frequency, mode intensity, as well as 
displacement pattern for the 105 infrared active modes. 
The calculations, which utilized the structure and valence 
as reported by Becker et al.^ were based on a real-space 
summation of screened coulomb interactions involving a 



spherical cut-off boundary^ The following two sections 
describe the method used to compare the observed and 
calculated spectrum for the A u and B u modes respec- 
tively. 
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B. 52B U modes 



modes at 20 K are found in Table Hill 



All 52B U modes predicted lie in the ac plane and 
should possess a finite intensity along both measured po- 
larizations in this plane unless the direction of the in- 
duced dipolc moment for a particular vibration is par- 
allel to either Sa or Be- (Small angle deviations from 
the parallel case may be hard to detect.) For all other 
modes one must sort out which excitations along e-A and 
ec corresponding to the same physical vibration so that 
they are not counted twice. To count the modes, we 
used the Drude-Lorentz parameters for each phonon at 
20 K as well as the resonant frequency and mode inten- 
sity from our lattice dynamical calculation. With knowl- 
edge of the optical conductivity spectra along ca and ec 
at 20 K, modes were grouped together if they were ob- 
served in the same small frequency interval (on the order 
of a typical linewidth) of one another in the two spec- 
tra. Next, we compared the experimental linewidths of 
the grouped modes to ensure that they would overlap if 
plotted side-by-side. As two final criteria, we compared 
both the ratio of experimental oscillator strengths with 
the computed oscillator strengths, and also compared the 
displacement patterns generated by our calculation to 
make sure they were in good agreement as to which exci- 
tations corresponded to the same physical vibrations. All 
52B U modes predicted in the ac plane were observed in 
the experimental spectrum. Lorentz oscillator parame- 
ters for all modes as well as modes present along both &a 
and ec polarizations can be found in Table U and Table 

m 



C. 53A U modes 

Counting the experimentally observed A u modes is 
straightforward; we observe 43 of the 53 predicted modes 
along this direction. We believe that 10 modes were not 
resolved as a result of weaker modes being buried within 
stronger modes and merging together in the spectrum. 
In the optical conductivity, all single phonon excitations 
should be represented by symmetric lorentzian oscilla- 
tors. Figure [7] left panel depicts a slightly asymmetric 
lorentzian oscillator with a resonance of 85 cm -1 . Our 
attempt to fit the resonance at 85 cm -1 with a sym- 
metric lorentzian oscillator yields unaccounted for spec- 
tral weight on the high frequency side of the resonance. 
A symmetric lorentzian oscillator centered at 79 cm -1 
and corresponding fit are also included as a reference. 
(Note that the additional spectral weight between the 
resonances at 75 cm -1 and 85 cm -1 is attributed to 
absorption between the bands and therefore not to be 
mistaken for a lattice vibration). We suspect that this 
asymmetric shoulder is due to a weak buried mode that 
is not fully resolved. To further support our interpre- 
tation of buried modes, the right panel of Fig. [7] depicts 
the appearance of a once buried mode (^276 cm -1 ) upon 
cooling. The Lorentz oscillator parameters of the 52_B U 
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FIG. 7. (Color online) The left panel (a) shows unaccounted 
for spectral weight on the high frequency side of the asym- 
metric oscillator at 85 cm -1 along with a fit to a symmetric 
oscillator (79 cm - ). We suspect the asymmetric shoulder is 
the result of a buried mode. The right panel (b) depicts the 
appearance of a once buried mode (~276 cm - ) upon cooling. 



D. Assignment of Modes 

The assignment of a particular phonon, observed in 
the experimental spectrum, to its precise atomic displace- 
ment pattern cannot be made without sacrificing objec- 
tive certainty because of the density of the phonon spec- 
trum. Assignments to groups of modes, however, can 
be made through the displacement patterns generated 
by our lattice dynamical calculations. Using a simplified 
model of a phonon as a driven oscillator with resonant fre- 
quency inversely proportional to mass, we proceed to list 
the constituent atoms in the structure from heaviest to 
lightest: Te, Br, Fe, and O. Normally, the low frequency 
modes involve the heaviest atoms, but in the FeTe2C>5Br 
structure the Br atoms have a weak interaction with the 
other atoms and therefore are responsible for the lowest 
frequency phonons. At a certain frequency the motion 
of Br atoms will cease. Likewise, at a higher frequency 
the vibration of atoms involving Te will cease, leaving 
the motion of Fe and O to be expected in the highest fre- 
quency phonons. Accordingly, the phonon spectrum can 
be divided into three clusters that is reflected in Fig. [2] 
(namely, below 180 cm -1 , between 180 cm -1 and 530 
cm -1 , and above 530 cm -1 ). The clustering is also ob- 
served in the calculated spectrum via the displacement 
patterns produced, in terms of first the Br and subse- 
quently the Te ceasing to vibrate. Within each cluster 
trends are to be noted on the following principle: atoms 
are listed corresponding to their net atomic displacement. 
This approach, in turn, has allowed us to identify an- 
other characteristic of atomic displacement within the 
three clusters. The atom of focus in a cluster (e.g., Br 
in the first cluster) has a net atomic displacement ever 
diminishing. Cessation does not arrive abruptly. The 
three clusters are profiled in each of Table IIII1 Table |H 
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and Table HI 

Apropos the third cluster that focuses on Fe and O, an 
additional resource is provided by the familiar orthofer- 
rite family of compounds because the FeOe octrahedral 
building block that is present in FeTe2C>5Br is itself a 
defining unit of the orthoferrite family. Given the abun- 
dant literature on vibrational assignments made to the 
FeC>6 ochrahedral unit ) 22 1 23 we are consequently able to 
map these well recognized assignments onto our experi- 
mental findings, thus corroborating both perspectives. 

V. CONCLUSIONS 

Far-infrared measurements on multiferroic single- 
crystal FeTe2C>5Br reveal no drastic anomalies in the ex- 
citation spectrum below the multiferroic ordering tem- 
perature. Nonetheless, a thorough inspection of the 
phonon spectra has lead to the identification of all 52 
predicted modes in the degenerate ac plane and 43 of 
the 53 predicted modes along the unique b axis. The ab- 
sence of 10 modes along the b axis is presumably due to 
the merging of weaker modes near stronger modes in the 
dense excitation spectrum. 

The motions of individual atoms have been assigned 
to groups of phonons via the displacement patterns pro- 
duced through our lattice dynamical calculation. The as- 
signments divide the overall phonon spectrum into three 
clusters, which are characterized by an atom that is com- 
mon to every vibration in that cluster. We observed that 
as frequency increases, vibrations involving the motion 
of Br atoms will progressively cease. Likewise, at higher 
frequencies the vibrations involving Te atoms will cease 
in the same manner, leaving the motion of Fe and O to 
be expected in the highest frequency phonons. Moreover, 
the particular atoms involved in a vibration are listed in 
descending order of net atomic displacement, which in 
turn makes clear a trend within each cluster (namely, 
the atom of focus in a particular cluster has a net atomic 
displacement ever diminishing). 



It remains to compare our results to a recent infrared 
study by Pfuner et ah 11 on the similar FeTe2C>5Cl com- 
pound where an equal number of phonons were predicted. 
Our analysis of the phonon spectrum agrees with the in- 
frared study of Pfuner et al. in that we also believe modes 
resonating at equivalent energy scales can overlap, pre- 
venting them from being resolved separately, but we dif- 
fer on other issues. Pfuner et al. report a combined 38 
modes along the e~B and ec directions. We observe more 
modes along the unique b axis alone, as well as modes 
below 50 cm -1 , whereas the study by Pfuner et al. does 
not report measurements below this frequency 

We also observe drastically different and more strongly 
anisotropic high frequency optical properties. Pfuner 
et al. report an absorption edge around 4,000 cm -1 
in FeTe20sCl, where as we observe low absorption 
accompanied by substantial transmission through the 
FeTe205Br crystal in this region. In speculating about 
the differences in high frequency optical properties of 
the two similar compounds, we point out that spurious 
changes of slope of the reflectance throughout the mid- 
infrared region, such as those resulting from the contri- 
bution of back surface reflection (black line in Fig. [5]), 
can drastically change the Kramers-Kronig-derived opti- 
cal properties. 

Transmission measurements as well as reflectance 
along &a were not reported in the study by Pfuner et 
al.. Similar to the FeTe2C>5Br system reported here, the 
FeTe205Cl system also showed no response to magnetic 
fields above or below the magnetic ordering temperature. 
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